Evolution of a sandpile in a thick flow regime 
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We solve a one-dimensional sandpile problem analytically in a thick flow regime when the pile 
evolution may be described by a set of linear equations. We demonstrate that, if an income flow is 
constant, a space periodicity takes place while the sandpile evolves even for a pile of only one type 
of particles. Hence, grains are piling layer by layer. The thickness of the layers is proportional to 
the input flow of particles ro and coincides with the thickness of stratifled layers in a two-component 
sandpile problem which were observed recently. We find that the surface angle 6 of the pile reaches 
its final critical value {6f) only at long times after a complicated relaxation process. The deviation 
{9f — 9) behaves asymptotically as {t/ro)~^^'^. It appears that the pile evolution depends on initial 
conditions. We consider two cases: (i) grains are absent at the initial moment, and (ii) there is 
already a pile with a critical slope initially. Although at long times the behavior appears to be 
similar in both cases, some differences are observed for the different initial conditions are observed. 
We show that the periodicity disappears if the input flow increases with time. 

PACS numbers: 83.10.Hh, 83.70.Fn, SS.lOPp, 46.10.-hz 
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I. INTRODUCTION 

Granular flows have attracted increasing interest over 
the last years, and constitute now a very active research 
field with apparent technological applications . Rea- 
sons for such great attractiveness of granular media or of 
the so called soft matter are clear: they present not only 
unusual properties which let them look like solids in some 
cases and liquids in other ones, but they also display new 
phenomena unknown both to solids and liquids. Such 
a new and intriguing phenomenon is recently observed 
spontaneous stratification of granular mixtures [p^-p^. 
When a granular mixture is poured into a Hele-Shaw cell 
(i.e. in a quasi- two dimensional silo) under some con- 
ditions, layers of different grains appear to be arranged 
periodically. 

Several studies were carried out to explain the strik- 
ing stratification |10-16| although the problem seems to 
be very complicated for an analytical treatment. In our 
communication we put forward the following ques- 
tion: is it possible to find some precursor of stratification 
in a much simpler situation when only one type of grains 
is poured? We shall present a full answer below for one of 
the limit regimes of sandpile spilling using an approach 
most convenient for analytical consideration. 

Crucial progress in understanding of granular flow na- 
ture was made some years ago when Bouchaud, Gates, 
Prakash, and Edwards |jl^ (see also the papers of Mehta 
and Barker |^9| and of Mehta, Needs, and Dattagupta 
|pO| ) introduced a clear phenomenological theory that de- 
scribes a surface flow of granular materials. The main 
idea of the approach is the following. The grains are di- 
vided to two parts - static grains and rolling ones. The 
conversion of static grains to rolling ones occurs depend- 
ing on the relation between the local slope of a pile and 
the repose angle of the material. In addition, the con- 



tinuity equation for the total amount of grains holds. 
Thus, the problem can be reduced to two coupled par- 
tial differential equations for the local amounts of static 
and rolling grains. One should note that similar ideas 
were also applied to a self-organized criticality problem 

Nevertheless, even this simple approach is still compli- 
cated because of nonlinearity of the equations. Recently, 
Boutreux, Raphael and P.-G. de Gennes |Q proposed 
a phenomenological description of some special case of 
granular flows - a so called thick flow regime - that pro- 
vides a unique possibility for an analytical treatment: 
in this case the coupled equations are linear and one of 
them is completely decoupled from the other one. In pa- 
pers ||2^,|2^ this approach was used to study some gran- 
ular flow configurations. Relaxation of sand from sev- 
eral most simple states was considered, generalization to 
many-component flows were proposed |26|-|2q| , and new 
improvements of the approach were made |29f| to make it 
closer to reality, though the main classical sandpile prob- 
lem still remained unsolved even for the simplest version 
of such a description. 

We shall use the proposed phenomenological equations 
to answer positively the above stated question and to de- 
scribe a total evolution of the sandpile, that will turn to 
be surprisingly complex. We will show that the slope of 
the pile approaches its critical value only at long times 
after a complicated discontinuous relaxation process [0 . 
In fact, we shall reconsider a classical sandpile problem 
using, maybe, the simplest possible idea without appeal- 
ing to more refined approaches like, for instance, self- 
organized criticality |30 3^,^. 

To start with, let us write out the phenomenologi- 
cal equations we shall use. Boutreux, Raphael and P.- 
G. de Gennes Q describe phenomenologically the one- 
dimensional thick granular flows (the flow thickness sup- 
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posed to be much higher than the grain size) by the fol- 
lowing equations: 



dr 
dt 



dr 
dx 
dh 
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dh 

dx 
dh 
dx 



(1) 



Here h[x, t) is a profile of the static part of the material; 
r[x, t) is the width of a moving granular layer; we assume 
that the flow is from right to left and all rolling grains 
are supposed to move with an equal velocity v. r{x,t) 
and h(x, t) describe completely the evolution of a pile. 
The sum of Eqs. (|l]), d{r + h)/dt - vdr/dx = 0, has the 
form of the continuity equation for grains, vr(x, t) is a lo- 
cal flow. The right hand parts describe the conversion of 
the static grains to rolling ones and vice versa depending 
on a relation between a local slope dh/dx and 9f. The 
meaning of w„ is the velocity of the uphill fronts as we 
shall see later. Usually Vu > v 14,2^]. 6*/ is a critical an- 
gle, or a so called repose angle, that is the angle to which 
the sandpile will evolve. The deviations of the local slope 
from Of are supposed to be small. Eq. (|^) is valid for 
r{x,t) > 0. If r{x,t) = 0, then the right hand parts of 
Eqs. (|l|) are supposed to be equal to zero. Therefore, 
the function Vu(r), which equals const -r in the thin flow 
regime [ p^ , now is a constant Vu for r > and is zero at 
r = 0. 

The physical reasons to introduce equations (|l]) in such 
a form are described in ||2j,^. One can see that, at a 
given point, static grains converse to rolling ones only 
if the local slope of the pile is higher than the critical 
angle, and rolling grains converse to static ones, if the 
local slope is lower than the critical angle. These linear 
equations are much simpler for an analytical treatment 
than the previously proposed nonlinear equations for a 
thin flow regime pj| , p6| -[2^ in which the characteristic 
velocity w„ in Eq. Im) is replaced by r{x, t) divided by a 
constant with the dimensionality of time. 

In fact, Eqs. (g) represent, maybe, a minimal model 
for the description of granular flows. The equations are 
really very simple: the second equation is independent of 
the first one. The general solution of Eqs. (|]) may be 
written immediately in the form p3 : 



r(x, t) = u{x + vt) 



V + v, 



-W{x - Vut) - VuOft, 



h{x, t) — w{x — Vut) + VuOft, r > 0; 



(2) 



h{x,t) — const. 



0, 



where u{x) and w{x) are arbitrary functions. We shall 
use it to describe the sandpile evolution phenomenologi- 
cally neglecting, as usual 1 24 2^ , possible near-front de- 
viations from the thick flow regime. In fact, in paper 
p5[ , moving fronts of granular flows were studied in in- 
termediate situations between the thick flow regime and 
the thin flow regime. The results j25j show that such a 
neglect is possible. 



It is convenient to use also the relation describing con- 
servation of number of grains that follows from the con- 
tinuity equation: 



dx[h(x, t) + r(x, t)] 



vrot+ / dx[h{x,0) +r{x,0)] 



(3) 



It follows directly from Eq. (|^) that the solutions con- 
sist of the part w{x — Vut) moving with the velocity u„ to 
the right, the part u{x + vt) moving with the velocity v 
to the left, and the growing (for h) or decreasing (for r) 
homogeneous background. Speciflc initial and boundary 
conditions, which we shall use in Sees. II and III, will 
result in u{z) and w(z) consisting only of linear parts. 
Therefore, the solutions will have breaks between linear 
parts. Some of them will move to the right with the 
velocity Vu and others - to the left with the velocity v. 
Coordinates of left fronts of the rolling and static grain 
distributions coincide and, because of the decreasing lin- 
early (for r) or increasing linearly (for h) homogeneous 
background, may move with velocities lower then v if the 
front of r(x) is not of jump-like form. One may check, 
using e.g. Eq. (^), that a jump- like front of the rolling 
grain distribution r(x,t) may move only with velocity v 
to the left. As it follows from Eq. (^), linear parts of 
r{x) and h(x) can move only parallel to themselves or 
are motionless. Thus, local slopes of the distributions 
can change only discontinuously if an income flow is con- 
stant in time. We shall show the picture of the evolution 
more clearly for particular examples in Sees. II and III. 

Thus, to describe the pile evolution we only have to 
solve linear equations with initial and boundary condi- 
tions. The boundary conditions are deflned from the 
condition on an input flow (we set it constant usually) 
and from consideration of a moving front of the pile - a 
moving boundary. In fact, that is the most difficult point 
in solving problems of this kind. 

There are several possibilities to choose initial condi- 
tions. In sections II and III we consider two most in- 
teresting and natural cases: (i) grains are absent at the 
initial moment, and (ii) there is already a pile with a crit- 
ical slope initially. We shall show that the pile evolution 
is a complicated process in both cases. The slope of a pile 
approaches its critical value only at infinite times after 
a long relaxation process during which areas with differ- 
ent slopes are present. At long times, the pile evolution 
looks very similar for both initial conditions, although, 
some distinction can be found even at the infinity. Thus, 
we observe a long memory of initial conditions. 

The results obtained for initial conditions (i) and (ii) 
are presented in Figs, p] and ^ 
bution profiles, in FigsT^ and 
of fronts and breaks of the grain distributions on the i, x- 
plane, and in Figs. ^ and ^- the dependencies of the pile 
slopes on time. 



- movies of r and h distri- 
i - diagrams of trajectories 
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In fact, all these figures depict a long relaxation pro- 
cess to the critical state (i.e. to the pile with the critical 
angle) by different means. For example, from Fig. 3 one 
sees that the relaxation is discontinuous - local deriva- 
tives of the slope over time are zero or infinity for all 
times. 

We shall demonstrate the space periodicity appearing 
while a pile evolves: it turns out that grains are piling 
layer by layer (see Fig. ^ and ^) and the thickness of 
the layers just coincides with the thickness of the strati- 
fied layers observed in two-component sand-piles Jll| , p^ ! 
Thus, we give positive answer to the question stated 
above and find a precursor of an intriguing stratifica- 
tion phenomena already in an one-component sandpile 
problem. 

One may wonder now if there is any possibility for a 
sandpile to evolve without a space periodicity at all. We 
demonstrate that, of course, such evolution is possible in 
a thick flow regime, for example, if an input flow is a lin- 
ear function of time. We show that the front moves with a 
constant velocity all the time, and linear profiles of mov- 
ing and static grains expand without any surprises in that 
situation (see Fig. The relation between amounts of 
static and rolling grains in a pile depends on the rate of 
the input flow increase. 



II. SANDPILE EVOLUTION STARTING FROM 
EMPTY STATE 

A. Constant input flow 

The case we will consider in this section corresponds 
to the situation when there are no particles in the initial 
state. 

Let there be a wall at a; = 0, and grains be poured 
permanently beginning from the moment t = at this 
point, so r{x = 0,t) — ro is a boundary condition (rp is 
the thickness of the input flow). Let, first, ro be constant 
in time. The sandpile is supposed to expand to the left, 
i.e. to a; = — oo. There are no particles at the initial 
stage so r{x < 0, t = 0) = 0, h{x < 0, f = 0) = are the 
initial conditions. Inserting general solutions from Eqs. 
(^ and the above initial and boundary conditions one 
may obtain the functions u{z) and w{z). 

To clear up the picture and to simplify the calcula- 
tions we start from the following consideration of the be- 
gininning stage. (Below, we shall demonstrate consistent 
complete calculations in the more frequent case, which 
is realized at more long times, of a r[x) front linearly 
decreasing to zero.) 

One may see that because of the initial and boundary 
conditions for r{x,t) and of the form of general solution 
for it Eq. (0), r{x,t) should first has a jump front; the 
part of r{x,t) adjacent to a; = has to be independent of 
t. The last statement follows immediately from the form 
of the general solutions Eq. (||) and from the fixed right 



boundary condition r{x = 0,t) — tq. The first state- 
ment follows from the last one - the jump front of r{x, t) 
is the only possibility for the first stage of the process. 
But then, as we have seen in Sec. I, it has to move with 
velocity v. 

The distribution of static grains h(x, t) which appears 
from the precipitation of rolling ones has a front that goes 
to zero without an abrupt leap but linearly in the case 
under consideration, since there are no fixed boundary 
conditions for h in the right. Then one can easily imagine 
that at first the solutions have to look as shown in Fig. 
1(a) (as we have noticed, nothing but linear functions 
may appear from the considered initial and boundary 
conditions). In principle, one may immediately insert r 
and h in this form ( r{x, t) = rg — ax, h{x, t) = h{x + vt), 
where a and h are constant coefficients) into Eq. (|l|) 
and obtain the following answer in the time interval 
< t < ii = (u + Vu)ro/{vVuOf ): 

r(x,t) = I j'o H Ofx] e{x + vt) 

V V + Vu ) 

h(x,t) = 9f(x + vt) (4) 

V + Vu 

Q{x) is the Heaviside function (we do not write the mul- 
tiplier Q{x + vt) in the right hand side part of the second 
equation of Eqs. (|^), since, of course, r{x) and h{x) can 
not be negative). The static grains are precipitated from 
the rolling ones, so the front coordinates of r{x) and h{x) 
should coincide. The meaning of the time ti is clear from 
Fig. 1(a) - it is the very first instant at which the jump 
at the front appears to be zero. Note, that the slope of 
the static part dh/dx — VuOf/{v -I- t;„) is less than the 
critical slope, so the relaxation to the critical slope is non 
trivial. 

Nevertheless, to be sure, we prefer to substitute the 
general solutions Eq. (^) into the following set of initial 
and boundary conditions 

r{x,t = 0) = , 
h{x,t = 0) = , 
r{x = 0,0 < t < ti) = ro , 
h{x = -vt,0 <t <ti) = . (5) 

Using previous considerations to simplify the calcu- 
lations, we assume in Eq. (^ the front velocity to be 
equal to V from the very beginning, so a fifth condition 
for r{x, t) at the front point is not necessary. In princi- 
ple, the time ti may be obtained if one demands that all 
parts of functions u{z) and w{z) have to be connected 
together continuously, but we prefer to use the already 
known expression for ti and be sure that calculations are 
self-consistent (i.e. all parts of functions u{z) and w{z) 
are connected together continuously) only afterwards. 

Substituting the general solutions Eq. (^ into Eq. (|^) 
we get 



3 



VVu Of 



< z < = 



V + Vn 



fZ , 



U I < Z < — I = To H J ; ^T^OfZ. (6) 



{V + Vuf 



After substitution Eq. (||) into Eq. (H) we again obtain 
our solution Eq. (^) of Eq. (|l|) with a front moving to 
the left with the velocity v in the time interval Q <t <ti 
[see Fig. 1(a)]. 

As a result, at the time ti one gets 



r{x,ti) ro H -^9fX, 



h{x, ti 



(7a) 



(7b) 



for -[iv + Vu)/ivv^)]ro/ef <x <0 [Fig. 1(b)]. These 
equations are used as initial conditions to find the solu- 
tions in the next time interval ti < t < t^, at which the 
front of r{x,t) will be jumpless. 

Times ^2 and ^3 appear naturally from the solution 
(see Fig. 1) but we prefer to write out their expres- 
sions immediately: t2 = [{v + /vv'^rn/Of and = 
[{v + Vu){v -\- 3w,j)/OT^]ro/6'y. As one may understand 
from the figure, the meaning of the times is the follow- 
ing: t2 is the time at which breaks of r(x,t) and h(x,t) 
moving from the front to the right with the velocity w„ 
will approach the point a; = 0; and ^3 is the time at which 
the break of r(x, t) moving from a; = will overtake the 
front. 

To simplify our calculations, we again use the expres- 
sion for ^3 and check the correctness of the choice of 
at the very end of the calculations. 

Now we show how to treat moving boundary condi- 
tions. We have to add to the initial conditions Eq. (0) 
the following boundary conditions: 



r{x = 0,t > 0) = To, 



(8a) 



V + Vu rp 

VVu Of 



Vd{t~h),ti <t<t3] = 0, (8b) 



h(x^ -!i±^l£ _ v^{t - h), ii < i < i3 ) - 0. (8c) 
V Of J 

Eqs. (|8b|) and (|^ are the conditions for the left front 
of the pile that is supposed to move with a yet unknown 
velocity Vd- All we have to do is (i) to insert the gen- 
eral solutions Eq. (^) into Eqs. (|^) and (||); (ii) to find 
Vd, u{z), and w{z); and (iii) to check the self-consistency 
of our choice of ^3 which was made beforehand, in fact, 
from physical reasons. Note that the beforehand choice 
of ^3 was used only to make the following intermediate 
formulas more compact. 

Substitution of Eq. (|) into Eq. (|^, (0), (|^, and 
(|^) gives 



w{ ^{vdv + 2vuV + 2vdVu + ivl)— < z < 

wi Of 

{v -\- Vuf rp {v + Vu){v - Vd) 



. , , _^ , rp + ^—OfZ, (9) 

VVu Of v[Vu+Vd) Vu+Vd 



{v + VuY rp v + Vurp\ Vu „ 

W I — < Z < — I = tifZ, 



VVu tif 



V Uf 



(10) 



M < Z < 



v_+Vu_rp_ 

Vu Of 



Vu{v + 2Vu) , . 

= '^°+ iv + vuy ^''^ 



and 



u < z < 



{v - Vd){v + Vu){v + 2vu) rp 



vvt 



7 w \0f^ 

[V - Vd)(V + Vu) 



(12) 



correspondingly. The functions u{z) from Eqs. (|ll|) and 
( p^ ) have to coincide since they are defined at (at least) 
overlapping intervals of z (we shall see below that, in 
fact, these intervals coincide). Thus, equating Eqs. (|Tl| ) 
and ( p^ we obtain the answer for Vd in the case under 
consideration: 



Vd = 



{v + 2vu) 



(13) 



One sees immediately that Vd < v,Vu. Inserting this ex- 
pression into Eq. ( [l2|) we find that the intervals of the 
variable z for u(z) in Eqs. ( [TTl ) and ( p^ are the same. 
Inserting Eq. ( [l3| ) into Eq. (||) we obtain 



^{v + Vu)'^rp {v + Vu)'^rp 
lu —3 — < z < — 



VVu 



V + Vu , {v + 2vu). 

rn H rOfZ. 

2vu 2{v + Vu) ^ 



(14) 



At last, substitution of Eq. (g) into Eq. ( |8a| ) gives, ac- 
counting for already known answers for w{z) Eqs. ( p^ 
and (|lj) in different intervals of z, the following expres- 
sions for u{z): 



\v + Vu)rp {v + Vu) rp 
u I — < z < T — 

Vu{v + 2Vu) 



ro + 



and 



[v + Vu)'^ rp 



< z < 3 



-.ro 



{V + U„)2 

{v + VuY rp 
vl Of 

Vu[2v + iVu) 
2{v + VuY 



(15) 



(16) 
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Combining Eqs. ( |ll[ ) and ( |l5| ) we get u{z) and One may see easily that they are bounded 

/ ( , \2 \ ( A. o \ continuously, so our choice of ^3 was correct. Inserting 

u(0<z< l!i±i^^ = ro + ')^''^ f OfZ. (17) Eqs. @, (0), ® , and Q into Eq. (|) we obtain a 



+ fuU solution of the problem for < i < ig: 

Eqs. (|l|), ®, and (|l|) present a fuU answer for 



r(vu(t -t2) <x <0,t>ti) = ro + ^" gfx, (18a) 



+ (18b) 



r(-w(i - ta) < a; < 0, t < tg) = ro + ^/x, (18c) 

h{vu{t - 12), ~vt<x<Q)^ ^" 9f{x + vt), (18d) 

[V -\- Vu ) 

-i < x < -Vu{t - ta), 0, i < <3 = ro + —, — ■ r Of [xA ■ (18e) 



Vu{v + 2vu)0j V + 2vu ' ' / 2vu 2{v + Vu) \ v + 2vu 



These solutions are shown in Fig. l(b-e) [see also Fig. 

the region ti < t < ^3 in which trajectories of the front 
- a lower segmented line - and the b reak s of r( x, t) and 
h{x,t) are depicted]. The solutions (18a) and ( [l8d| ) are 
defined in region 1 of Fig. the solution (|18b|) is defined 
in region 2 of the figure, th e so lution (|l8q ) is defined in 
region 3, and the solution (18e) for h{x,t) is defined in 
both regions 2 and 3. 

The evolution looks as the following. At the instant 
ti new linear parts of r and h appear at the front point 
X = —[(v + Vu)/vu\ro/9f. The new front moves to the left 
with the velocity Vd, breaks of r(x, t) and h{x, t) move to 
the right with the velocity - a solid line x = Vu{t — t2) 
in Fig. 2. The slope of this part of the static grain dis- 
tribution h{x, t) is higher than it was at the first stage 
but lower than the critical slope. The breaks approach 
the X = Q point at t = ia- Then a new time-independent 
linear part of r(a::) appears close to the wall. The left part 
proceeds to move to the left, so we see a break moving 
downhill - a dashed line x — —v{t — ^2) m Fig. 2. Its 
velocity equals v. Note that, unlike r{x,t), h{x,t) has 
no breaks for t2 < t < t^. At the instant ^3 the break 
of r{x, t) overtakes the front. That is possible, since the 
velocity of the front is lower than the downhill velocity 
of the break, Vd < v. 

Therefore we confirm the general picture of the evolu- 
tion predicted in Sec. I: there exist breaks of both r{x, t) 
and h{x, t) moving uphill with the velocity u„, and only 
the break of h{x,t) is moving downhill with the veloc- 
ity V. When the breaks approach a wall or the break of 



h{x,t) overtakes the front, new linear parts of the solu- 
tions appear. 

From our solutions Eq. (p^), we obtain new initial 
conditions for r{x,t) and h{x,t) at the moment t = tg. 
Then we may repeat the described procedure for the next 
time interval ^3 < i < is, etc. However, now, when we 
understand the structure of the solution, one may again 
simplify the calculations. Instead of calculating the func- 
tions u{z) and w{z) one may proceed directly with r{x^ t) 
and h{x,t). Areas of different linear parts of solutions 
are triangular regions in Fig. ^ separated one from each 
other by solid or dashed lines - trajectories of the breaks. 
Let us suppose that the solution is known in the region 
(0, 0) — (xi, ii) — (0, ta), i-e. in region 1 in Fig. ^. One may 
sew easily an unknown linear solution in the adjusted tri- 
angular region 2 together with the known one along the 
line x = —v{t — ia)- It is easy to find coefficients of linear 
terms and, therefore, the front velocity. From this values 
one finds the shape {xi,ti) — (0, ^2) — {x2, t^) of region 2 
in which the new part of the solution is defined (see Fig. 
H). Then we repeat the described procedure for the next 
triangle, etc. 

After these simple but rather tedious calculations we 
obtain the total solution consisting of linear parts, a 
structure of which one can see from Figs. 1 and ^ and 
following equations. The lowest segmented line in Fig. 
^ shows the dependence on time of the front position. 
Coordinates of the segments are: 



x(t 



2m- 1 



<t<t 



2m+l 



m(m + 1) {v + Vu)'^ ro 

2 VuVm Of 
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VUu 



t, m = 0,1,2, 



(19) 



(for m = the time interval is < t < ti) with 
Vm = Tnv + (m + l)vu- The particular times shown in 
Fig. 2 are 

TO (w + Vu)[[m - l)v + (to + l)Vu] To 



t2 



m — 1 



m{m + 1) {v + Vu)^ To 



(20) 



TO = 1,2, 



The front coordinates corresponding to times ^2™-! at 
which the front velocity changes its value are 



X[t2m-l) = -m — ,TO=1,2... (21) 



Hence, these points are arranged periodically. Two other 
types of lines are shown in Fig. 2: solid lines x — 
Vuit — t2m) and dashed lines x = —v{t — t2m), to = 1,2,... 
The lines of the first type depict the uphill movement of 
the breaks of both profiles r{x) and h{x) [Fig. 1(c)]. The 
lines of the second type show the downhill movement 
of the break of the profile r{x) with the velocity v [see 
Fig.l(e)]. 

Solutions for all regions of Fig. 2, connected at these 
lines, look as: 



r {-v{t - t2m-2),Vu{t - t2m) < ^ < 0) = Tq + 



m{m + 1) [v + Vu)^ ro 



VuVm Of V, 



——t < X < -V[t - t2m.),Vu{t - t2m) = TT 



to(to + 1) (w + Vu) 



VuVm „ ( VV. 
H t I X -\ 

TO = 1,2, 



Ofx, 



eJx+^^t), (22) 



h{Vu{t - t2,n+2) <X < Vu{t - t2„),0) = — -Tq 



(m + l){v + Vu) 



Here, we set = 0. Inequalities in the right hand parts 
of Eqs. ( p^ ) and (|2^) define the areas of validity of the 
solutions (see Fig. 2). Eqs. ( p2| ) and ( [23| ) describe totally 
the sandpile evolution, see Fig. l(b-e). 

From Eq. ( p3| ) and Fig. 2, a space periodicity of the 
process is evident: general shapes of the profiles r(x) 
and h[x) are repeated each time the front moves by 
{v -\- Vu)'ro / {vu6 f) to the left. In fact, the pile is increased 
layer by layer, and the expression for the width of these 
layers is the same as the one for the width of stratified 
layers of different fractions in a two-component sandpile 
(that was obtained in papers [|l^,|lj]). One sees from Eq. 
(p3|) that heights of the pile at times t2m are also periodic 



h{x^O,t2m) = ni'^^^ro ,to = 0,1,2... (24) 



It follows from Eq. ( |23| ) that the slope, 9h = dh/dx, of 
the static part of the pile will approach its critical value 
only at infinite time: 



,TO = 1,2. 



')h{t2m-3 <t< t2m) — 1 ^ 



TO (w + Vu) 



(25) 



0, 1,2, . . . all the profile has the same slope. Thus, for 
long times t » [{v + Vu)'^ / {vVu)\rQ / 9 f the slope behaves 
as 



2v6ft 



(26) 



and Oh relaxes slowly to its final value 0/ by a power law. 

A slope of the whole pile - including both static and 
rolling parts, 6*;,+,. = d{h + r)/dx, depends on time in 
the following way: 



Oh+rit2,n-2 < t < t2m) = ( 1 + ^" " ^ ) Of 



W (w + Vu) 



r(i2m-l < t < t2 



m{v + Vu) ' m(rn + l){v + Vu)"^ ' 
m = 1,2... 



(27) 



(see Fig. ||) and approaches the final value f from above 
at long times. It behaves asymptotically as 



(for TO = 1, < t < i2) (see Fig. 3). At t2m-i < t < 
t2m,TO = 1,2,... there are two different slopes for two 
parts of the profile h{x). For t2m < t < t2,n+i,TO = 



[compare with Eq. (26)]. 

At long times, the coordinate of the front is x = 
y/2vrot/9f, and its velocity tends to — vtq / [29 ft) . Two 
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last relations are obvious. Indeed, at long times the slope 
of a pile is close to its critical value. The relative part 
of rolling grains is small. After substitution of a lin- 
ear function with the critical slope into the condition of 
conservation of grains Eq. (j^) we obtain last relations 
immediately. 



B. Time dependent input flow 

As we have noted in Sec. I, space periodicity of the 
pile evolution may disappear if an input flow increases 
with time. Let us consider the simplest example of a lin- 
ear time dependence: ro(t) = vqI (the constant vq has 
dimensionality of a velocity). Now there are not any 
jumps of profiles at the front point, since the growth of 
ro starts from zero and an initial state is empty, so the 
problem is simpler than that considered above. 

Instead of trying to solve the problem directly, one may 
proceed in the following way. First we suppose that the 
velocity of the front is constant. Let us call it Vd again. 
From a moving boundary condition h(x — —Vdt, t > 0) — 
at the front point accounting for Eq. ^ we obtain 
w{z). Then from another moving boundary condition 
r(x = —Vdt.t > 0) = and Eq.(||) one may get u{z). In- 
serting both answers into the condition r{x — Q,t) — v^t 
we (i) verify that the choice of a constant front velocity 
is right and (ii) find Vd as a function of vq. 

We write out the answers immediately and check them 
by substituting in Eq. (^. The profiles of moving and 
static grains look as 

r = —{x + Vdt) , 

Vd 

h ~ — 0f{x + Vdt) 

Vd + Vu 

(see Fig. 4), where the front velocity is 



approaches its highest possible value, and rolling grains 
have no time to converse into static ones. For a low c 



(29) 



1 



Vd 



2 1 + OfVu/vo 



{Vu - t')^ 



1 + ej- 



vo 



(30) 



In the figure, the evolution of the pile in the situation 
under consideration is depicted. One may see that the 
relation between amounts of rolling and static grains de- 
pends considerably on the rate of the input flow increase. 
As it follows from Eqs. (p9|) and (jsO), when vq » 



(31) 



Vd = V, 

r = — (X + vt) , 

V 

—^ef(x + vt) , 

V + Vu ^ 



and the amount of rolling grains is much higher than the 
amount of static grains [Fig. 4(a)]. The front velocity 




(32) 



so the front velocity tends to zero, the slope of the distri- 
bution of static grains is nearly critical, and the relative 
amount of rolling grains is small [Fig. 4(b)]. 

The amounts of rolling and static grains in a pile are 
equal for the following particular value of the rate of the 
input flow growth: 



VUu 



" v + 2vu ' 

(compare with Eq. ( ^ from Sec. II A) 
lar case, the front velocity is — v/2. 



(33) 

In this particu- 



III. SANDPILE EVOLUTION STARTING FROM 
CRITICAL STATE 

Now, let us switch on the input flow of grains when 
there is already a pile with the critical slope. Thus, the 
initial condition for h is h{x,t = 0) = 9f{x + d), where 
d is the initial horizontal size of a pile. All other initial 
and boundary conditions are the same as in Sec. II A 
(the input flow is time-independent). 

Now one can proceed with the calculations similar to 
those that were made in Sec. II to obtain the solutions 
describing the evolution of the pile. We omit the te- 
dious calculations and present immediately results for 
the case under consideration. Main answers are pre- 
sented schematically in Fig. 5 in which the distributions 
of rolling grains r{x) and the static ones h{x) are shown 
in several successive moments. (Note that scales change 
from one flgure to another.) Trajectories of the front and 
the breaks on the t, a;-plane are shown in Fig. 6. 

Let us comment the contents of Fig. 5, since it looks 
rather intricate. At the initial moment rolling grains are 
absent, and the angle of the profile of static grains is criti- 
cal. After we switch on an external fiow, a step of rolling 
grains starts to descend downhill the critical profile of 
static grains [Fig. 5(a)]. Rolling grains do not convert to 
the static ones. 

After the step approaches the last point of the pro- 
file of static grains [Fig. 5(b)] it proceeds to move to 
the left with the same velocity v, but the height of the 
jump starts to decrease with time [Fig. 5(c)]. A new lin- 
ear part of r{x) emerges and begins to move uphill with 
the velocity Vu- (Recall that all linear parts of profiles 
move as a whole or stay without movement. They can 
not change their slope with time.) A new more gently 
sloping part of h{x) simultaneously appears and begin to 
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move to the left with the velocity v, so a. break of h{x), 
coinciding with the break of r{x), moves uphill with the 
velocity Vu- 

After the decreasing jump at the front of r{x) disap- 
pears [Fig. 5(d)] the hnear part of r{x,t) proceeds to 
move uphill with the velocity [Fig. 5(e)]. The corre- 
sponding part of h{x) also proceeds to move uphill. Be- 
hind it, a new static part of the profile with the critical 
slope emerges, so, in fact, an inclined step moves uphill. 

When its first point touches a wall [Fig. 5(f)], a new 
static linear part of r(x) appears to satisfy the bound- 
ary condition r{x = 0,t) = ro [Fig. 5(g)]. An old part of 
r{x) proceeds to move uphill, so a new break of r{x), that 
moves downhill with the velocity v, emerges. Evolution 
of h{x,t) proceeds as before. At some instant [Fig. 5(h)] 
the last point of r{x) and the break of r{x) meet. New 
linear part of r{x) emerges and starts to move downhill 
[Fig. 5(i)]. The h{x) profile evolves in the same way as 
before. 

When all the h{x) profile appears to be critical [Fig. 
5(j)], the whole r(x) distribution is linear. This inclined 
profile starts to move down the critical one [Fig. 5(k)]. 
The only difference with the instant depicted in Fig. 6(a) 
is that the step of rolling grains is inclined now. After 
its first point approaches the last left point of the critical 
pile [Fig. 5(1)] new linear parts of r{x) and h{x) emerge 
which move to the left with some velocity Vd < v,Vu 
[Fig. 5(m)]. Two breaks of r(x) emerge which move one 
to each other. 

After the breaks meet [Fig. 5(n)], two new breaks of 
r(x) emerge which move away one from each other [Fig. 
5(0)] while h{x,t) proceeds the previous evolution. The 
moving downhill break of r(a:) overtakes the front at some 
moment [Fig. 5(p)]; after that an inclined step starts to 
climb uphill leaving after itself a new static critical slope 
part [Fig. 5(q)]. Afterwards the evolution proceeds in 
the same way. 

Thus, one sees that the evolution of a pile in this case 
looks more complicated than in Sec. II. Nevertheless, 
as it is evident from Fig. 6, in which all trajectories of 
the front and breaks of r(x, t) and h{x, t) are shown, a 
general similarity remains. 

Let us describe Fig. ^. The solid lines x — Vu{t 



^2m)' ^uit - t2m): 



0, 1, . . . show trajectories of the 



r{x,t) and h{x,t) breaks moving uphill (we introduce 
these notations for the instants of time to keep a tie with 
the corresponding notations in Sec. II). The dashed lines 
X = —v{t — t2m)i —v{t — t2m), TO = 0, 1, . . . and X = —vt 
depict trajectories of the r{x,t) break moving downhill. 
Comparing with Fig. 2 from the previous section one sees 
that the corresponding lines for the break trajectories are 
splitted now, and static segments of the front trajectory 
appear. The expressions for particular times shown in 
Fig. 6 are 



v + Vu,^ {m - l)m {v + VuY 

t2rn-2 = TO d H ^ 



v + Vu, mim + 1) {v + VuY 
TO a + 



2 vvl 
m = 1,2, , 



(34) 



- for the instants when the breaks of r{x,t) and h{x,t) 
reach a wall - and 

hrn-l = a + — 5— (to - 1)V + (to + 1 

vvu ^ vvi Of 

V™ TO + 1 W + W„ , 1 ^0 

Wu 2 vvf^ 0' 

TO = 0,1,, 



(35) 



- for the times at which the break of r(a;, t) overtakes the 
front. 

The trajectory of the front can be described by the 
following relation 



^(*2m+l <t < <2m+l) 



(to + 1) 



Vu Of ' 
X{t2m-1 <t< i2m+l) = 
to(to + 1) {v + W„)^ ro VVu 



VuVra Of V„ 
TO = 0, 1, 



(36) 



and a;(0 < t < t_i) = —d. Note that the second relation 
of Eq. ( p6| ) coincides with Eq. (p^). Thus, at the cor- 
responding instants, the front velocities are the same in 
both cases. The first relation of Eq. ( ^6|) again demon- 
strates the space periodicity of the process. 

Time intervals for which the front is motionless do not 
change with time: 



2m-l 



to 



V + Vu 



VUu 



d , TO = 1, 2, , 



(37) 



These intervals appear because of a new initial condition, 
so even at long times one may observe the influence of the 
initial condition. Nevertheless, the duration of the other 
time intervals, when the front moves, increases with time. 
Thus, the relative effect of the initial conditions decreases 
with time. 

We shall need only the following result for h: 

h{Vu{t - t2™+2) <X < Vu{t - t2m+2)) = 



Of \ x + d+{m + l)- 

\ Of 
h{Vu{t - <2m+2) <X <Vu{t- t2m)) 



m {v + Vu) 



-ro 



Vn 



(to + l){v + Vu) 



X+^-^t 
TO = 0,1,. 



(38) 



One sees that these relations are valid for different strips 

Vuit - t2m+2) < X <Vu{t - t2m+2) and Vu{t - ^3^+2) < 

X < Vu{t~-t2m) of the a;,t-plane. The first relation of Eq. 
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( p8| ) describes parts of the distribution of static grains 
with the critical slope, and the second relation of Eq. 
( ^ ) coincides with the corresponding answer Eq. ( |2^ ) 
for the first case. 

The slope angles of the profiles of static grains obtained 
from Eq. ( ^8| ) are shown in Fig. 7 which demonstrates 
most naturally the relaxation to the critical state. Unlike 
the previous section (see Fig. 2), in this case, segments 
with the critical angle are present. As can be seen from 
Fig. 7, different segments may overlap, and the profile 
may have two or even three parts with different slopes 
simultaneously (see also Fig. 6). Long time asymptotes 
of the pile slope are the same as in section II. 



eral picture of a sandpile evolution will hardly change, 
since the structure of the evolution equations is the same 
- in fact, the second equation of Eqs. (|l|) is independent 
of the first one. 

The following questions remain open. How crucially 
do our results depend on the used, maybe oversimplified, 
phcnomcnological approach? Are they really of general 
significance? Does the sandpile evolution in a thin fiow 
regime differ essentially from the described one? It is 
apparently impossible to answer these questions within 
the present approach, however the last results and 
psf on granular ffows in thin and intermediate regimes 
let one hope to solve these beautiful problems. 



IV. CONCLUSIONS 



In summary, in the case of the thick flow regime, we 
have shown that space periodicity takes place during a 
sandpile evolution even for the one-component pile. The 
pile spills during a repeating process if an income flow is 
constant: grains are piling layer by layer. The thickness 
of the layers coincides surprisingly with the thickness of 
stratified layers at the two-component sandpile problem 
0,0. Thus, in the one-component pile, we found a 
clear precursor of the stratification phenomena. 

We have found very reach behavior using the most sim- 
ple and clear approach |Q admiting an easy analytical 
treatment. The slope of the pile goes to its final critical 
value after complex relaxation which is a long discontin- 
uous process. Indeed, in the case of a constant income 
flow, the pile slope may change only discontinuously (see 
Figs. H and At long times the deviation of the slope 
from its critical value is proportional to ^/rojt. 

The pile evolution from different initial conditions 
looks very similar at long times but even at such times 
one can observe some difference (compare Figs. |^ and ^ 
for two most natural situations: the sandpile evolution 
starting from an empty state and the evolution starting 
from a critical state). 

The space periodicity disappears if an input flow grows 
linearly with time. In such a situation the sandpile front 
moves with a constant velocity all the time, and the pile 
evolves continuously without any peculiarities (see Fig. 
i- 

Formally speaking, we studied only the one- 
dimensional problem, but an ordinary three-dimensional 
sandpile with an axial symmetry (that means that sand is 
pouring on an infinite horizontal plane at a single point) 
may be described by the same equations as Eqs. (1) (the 
coordinate x means the distance from the center). Thus, 
our results also stay valid for such a pile if one applies 
the phenomenological description p^ , |2^. 

Aradian, Raphael, and de Gennes |29| have recently 
introduced a dependence of the downhill velocity v on r 
in the frames of the thick flow regime («« = const). The 
profile shapes in that case become nonlinear but a gen- 
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FIG. 1. The evolution of the profiles of rolling grains r{x) 
and static grains h(x) in the region — oo < a; < 0. In the 
initial state grains are absent, 

r(-(X) < I < 0,t = 0) = h{-oo < X < 0,t = 0) = 0. (a) 
< t < ti, the front moves with the velocity v. (b) t = ti 
[see Fig. 2 and Eq. (4)]. (c) ti < t < t2, the front moves with 
some velocity Vd < v < Vu, the breaks of the profiles move 
uphill with the velocity Vu- {d) t = t2, the breaks approach 
the wall at a; = 0. (e) t2 < t < ts, the front proceeds to 
move with the velocity Vd, the break of r{x) moves downhill 
with the velocity v. Note, that the right linear part of r{x) is 
always motionless. After the break of r{x) overtakes the front 
a,t t = ts, the general configurations (6) - (e) is repeated with 
a lower front velocity. 




FIG. 2. The areas of different solutions for the profiles of 
rolling and static grains for the pile evolution starting from 
an empty state [see Eqs. and (p3|)]. v/vu = 0.3. The 

lower segmented line shows the dependence of the front coor- 
dinate on time. The solid lines x = Vu{t — t2m),m = 1,2,... 
depict the uphill movement of the r(x) and h(x) breaks. The 
dashed lines x — —v{t — t2m) show the downhill movement 
of the break of the r{x) profile. The points Xm are arranged 
periodically. 
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FIG. 3. The dependence of the relative pile slope on 
time when the pile evolution starts from an empty state. 
v/vu = 0.3. Oh is the slope of a static part of the pile. 
9h+r is the slope of the whole pile consisting of static and 
rolling parts. The upper set of lines shows 9h+r- At infinity 
dh,Oh+r Sf- The separated lines for 6h+r are defined for 
< t < t2, ti < t < U, t2 < t < ti, etc. (see Fig. 2). The 
lines for Oh are defined for Q < t < t2, t\ < t < ti, ts < t < ta, 
etc. Sometimes (e.g., at ti < t < t2, < t < ti, etc. for Oh) 
the profiles have two parts with different slopes. 




X X 

FIG. 4. The evolution of a sandpile with a linearly growing 
input flow: ro = vot, vo = const. An initial state is empty. 
We show here both static and rolling grains in the same plots, 
so profiles for h{x) and h{x) + r{x) are presented. The fronts 
moves with a time- independent velocity Vd which is a function 
of Wo. (a) The rate vo is much higher than vq [see Eq. (^) and 
the text]. The front velocity is close to its maximal possible 
value V, rolling grains have no time to converse into static 
ones, so the amount of rolling grains is much greater than the 
amount of static grains. ( b) The rate vo is much lower than Vq . 
A relative amount of rolling grains is small, the front velocity 
tends to zero, and the slope of the distribution of static grains 
is nearly critical. 
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FIG. 5. The evolution of the profiles of rolling 

grains r{x) and static grains h{x) starting from the 
critical state (see the text for details). At the ini- 
tial moment r(— oo <a::<0,t = 0) = and 
h{-oo < X < 0,t = 0) = 9f{x + d). Plots (ffl)-(g) show 
the profiles at some successive instants. Afterwards the evo- 
lution proceeds in a similar way. The scale is changed from 
one plot to another. 
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FIG. 6. Trajectories of the front and the breaks of the 
r{x) and h{x) profiles on t, x-plane in the case of the pile 
evolution starting from a critical state (compare with Fig. 2). 
v/Vu = 0.3 . h{x,t = 0) = Of{x + d). The lower segmented 
line shows the dependence of the front coordinate on time. 
The solid lines x = Vu{t — tim), m = 1, 2, . . . depict the uphill 
movement of the r{x) and h{x) breaks. The dashed lines 
X — —v{t — t2m) show the downhill movement of the break 
of the r{x) profile. The coordinates of plateaus are arranged 
periodically. 
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FIG. 7. The dependence of the relative pile slope QhjQf 
on time when the pile evolution starts from the critical state 
(compare with Fig. 3). vjvu = 0.3. In some time intervals 
the profile of static grains has two or even three parts with 
different slopes. Segments with the slope equal exactly to the 
critical one are present even at long times. 
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